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We analyze the screening of an external Coulomb charge in gapless graphene cone, which is taken 
as a prototype of a topological defect. In the subcritical regime, the induced charge is calculated 
using both the Green’s function and the Friedel sum rule. The dependence of the polarization charge 
on the Coulomb strength obtained from the Green’s function clearly shows the effect of the conical 
defect and indicates that the critical charge itself depends on the sample topology. Similar analysis 
using the Friedel sum rule indicates that the two results agree for low values of the Coulomb charge 
but differ for the higher strengths, especially in the presence of the conical defect. For a given 
subcritical charge, the transport cross-section has a higher value in the presence of the conical 
defect. In the supercritical regime we show that the coefficient of the power law tail of polarization 
charge density can be expressed as a summation of functions which vary log periodically with the 
distance from the Coulomb impurity. The period of variation depends on the conical defect. In the 
presence of the conical defect, the Fano resonances begin to appear in the transport cross-section 
for a lower value of the Coulomb charge. For both sub and supercritical regime we derive the 
dependence of LDOS on the conical defect. The effects of generalized boundary condition on the 
physical observables are also discussed. 

PACS numbers: 81.05.ue, 03.65.-w 


I. INTRODUCTION 

Coulomb screening in graphene EMU provides a glimpse of the strong nonperturbative QED effects M and gives 
information about the measurable transport properties [Mi- The phenomenon_of the Coulomb screening can be 
analyzed in several ways. One method is to use the Green’s function technique [l5l - [20| . which was originally developed 
using the solutions of the corresponding eigenvalue equation fl5],|Ti|] . Subsequently a more powerful operator formalism 
for the Green’s function was introduced in QED [l7j-[l9| , which has been adapted to the case of planar graphene [20| . 
Another approach for studying Coulomb screening is provided by Friedel sum rule [2li [22| , which uses the spectral 
data in the scattering sector. In this approach, the induced charge is related to the scattering phase shift. This 
technique has also been applied to a variety of fermionic systems f23i - i2(I including planar graphene ana. 

If the external Coulomb charge in graphene exceeds a certain critical value, the system exhibits quantum instability. 
This is characterized by a rapid oscillatory behaviour of the wavefunction near the location of the charge 3, 0- 
In addition, quasi-bound states appear in the spectrum for gapless graphene 300 , for which there is recent 
experimental evidence [28[. These phenomena are manifestations of strong nonperturbative QED effects. 

Topological defects provide another source of nonperturbative quantum effects in graphene [29l - l38j . A graphene 
cone provides a simple prototype of a topological defect, which can be obtained by cutting a sector of planar graphene 
and gluing the two ed ges of the removed sector. This conical defect generates nontrivial holonomies of the quasiparticle 
wavefunction [29j - l3TI . I33l - l36i . HI, which can be modelled through the introduction of a fictitious magnetic flux tube 
with a suitable vector potential. The problem of the Coulomb charge in a graphene cone can therefore be related to 
an equivalent problem of the Coulomb charge in planar graphene with a suitably chosen magnetic flux tube. 

The main purpose of this paper is to show how the conical defect affects the Coulomb screening and related physical 
observables in gapless graphene. In the subcritical regime of the external charge, the screening is analyzed using both 
the Green’s function and the Friedel sum rule. We follow the operator technique [l8l - l20l | to calculate the Green’s 
function from which the induced charge density is obtained. Our calculation clearly shows the effect of the conical 
defect on Coulomb screening. In particular, a plot of the induced charge against the Coulomb strength indicates how 
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the conical defect modifies the value of the critical charge. We also calculate the induced charge using the Friedel sum 
rule [3j. For the case of planar graphene, these two approaches agree for low values of the external charge and starts to 
show divergence as the critical value of the charge is approached. In the presence of the conical defect, the difference 
shows up at much lower value of the external charge. This could be indicative of certain subtleties associated with the 
Friedel sum rule in the presence of singularities [391 ]. For both the Green’s function and Friedel sum rule approach, 
our calculations with the conical defect have the correct planar limit. For a given value of the subcritical charge, the 
transport cross-section has higher value in the presence of the topological defect. Similarly, the LDOS is also affected 
by the topological defect. 

In the supercritical regime, the_induced charge for planar graphene has already been obtained in the literature 
using both the Green’s function [40] and the Friedel sum rule jsj. Here we show how the topological defect affects 
the scattering phase shifts in the supercritical regime. This leads to the corresponding effect on the induced charge 
via the Friedel sum rule. The polarization charge density in the supercritical regime exhibits a power law tail. The 
coefficient of the power law tail can be written as a sum of functions which vary log periodically with the distance 
from the Coulomb impurity E3- We show that this period is affected by the presence of the conical defect. In the 
presence of the topological defect, the Fano resonances BE' appear for a lower value of the charge. This indicates 
that the conical topology affects the value of the critical charge. We also calculate the LDOS and find its dependence 
on the conical defect. These results could be of interest in the context of graphene based electronic devices [TTI. EUt 
In the supercritical regime also, our results have the correct planar limit. 

For certain parameter ranges in the subcritical regime, the gapless graphene system admits generalized boundary 
conditions |13l4lol |. When a Coulomb charge is embedded in a graphene sheet, it can induce various short range 
interactions. Effect of these interactions cannot be directly incorporated in the Dirac equation which is valid only 
in the long wavelength limit. The generalized boundary conditions encode an average effect of such short range 
interactions. The physical observables of the system depend on the choice of these generalized boundary conditions. 
Here we have shown how such generalized boundary conditions affect the LDOS and transport cross-section in the 
presence of the topological defect. 

This paper is organized as follows. In Section 2 we calculate the induced charge using the Green’s function in a 
subcritical Coulomb field in graphene cone. Next we analyze the subcritical region using the Friedel sum rule and 
show how the physical quantities depend explicitly on the sample topology. We discuss the effect of topological defect 
on the electrical conductivity of the system. We compare our results obtained using Green’s function technique and 
Friedel sum rule. In Section 3 the analysis of the corresponding spectrum is done in the supercritical region. This is 
followed by a discussion of the effect of generalized boundary conditions in the subcritical regime in Section 4. We 
conclude this paper with a summary in Section 5. 


II. SCREENING OF SUBCRITICAL COULOMB CHARGE IN GRAPHENE WITH CONICAL DEFECT 


A. Induced charge from Green’s function 


Our system consists of an external Coulomb charge impurity at the apex of a gapless graphene cone. In order to set 
up the formalism, we first consider planar graphene. For technical reasons [lsj, [Ej it is useful to introduce a mass A4 
for the quasiparticles in graphene. This serves as an infrared cutoff which is removed at the end of the calculations. 
The Dirac equation for a planar graphene in presence of a Coulomb charge is given by 


—i<Tid x — ioidy 


a 

- 1 - M(J 3 

r 


4 - = E 4 >, 


(1) 


where the wavefunction T has two components corresponding to the two different sublattices. The Pauli matrices 
ci, 2,3 act on the sublattice indices and a is the Coulomb potential strength. A positive value of a corresponds to a 
positive external charge. The radial coordinate in the x — y plane is denoted by r. The Dirac equation can be written 
in the form 


[Ho 


A4]T = 0, where Hq = a 3 E — (^—icrid x — ia- 2 d.. 


The Green’s function G{f,i T o\E) corresponding to Eq.([2j) satisfies the equation 



( 2 ) 


[Ho ~ M\G{r, ro; E) = 5{f- G 0 ) 
and the induced charge density pi(f) in graphene is given by 

J c 


( 3 ) 


pi(f) = — 4ie 


( 4 ) 
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where the factor 4 appears due to spin and valley degeneracies in graphene, e is the electronic charge, G(r, r; E ) 
denotes the Green’s function and C denotes the contour of integration. The Green’s function defined in Eq.® can 
be written as 

G(r, ro; E) = -i[H 0 + M][-ir{H 2 - M 2 }]' 1 ^S(r - r 0 )5(6 - 9 0 ), (5) 

where 6(r — ro) = -y= 6(r — ro)S(9 — 9q) in polar coordinates. Using the Laplace transformation 

pOO 

[~ir{H 2 - M 2 }]- 1 = / dse irs{ii *- M2 \ (6) 

Jo 

the Green’s function can be written as 

poo i—— 

G(r,f 0 -,E) = -i[H 0 + M} dse irs ^~ M2} J—S(r - r 0 )6(9 - 9 0 ). (7) 

Jo V r o 



FIG. 1: Formation of a graphene cone from a plane sheet of graphene by removing the sector POQ and identifying the edges 
OP and OQ. 

We now apply the above formalism to graphene with a conical defect. As shown in Fig.(l), such a defect is 
introduced in planar graphene by removing n sectors subtending an angle and then identifying the edges, where 
n can have values 1,2,3,4, 5. The value n = 0 denotes planar graphene. As a result of this identification, the frame 
{ e x , e y }, shown in Fig. Gib becomes discontinuous across the joining line. To take care of this discontinuity we choose 
a new set of rotated frames (30i. l3ll. l38i] e X ' = ee and e y > = —e r . The wavefunction associated with a quasiparticle in 
graphene picks up a holonomy when the quasiparticle encircles the conical defect. This is a topological effect, which 
can be modelled by considering a graphene plane with a fictitious flux tube perpendicular to the planar surface. The 
vector potential A associated with this fictitious flux tube is so chosen that when a quasiparticle encircles the flux 
tube, the wavefunction picks up exactly the same holonomy as in the case of the graphene cone. The resulting Dirac 
equation is given by 


Ot 

a.p' — eagAg + A4ct3- 

4/ = E \k where Ag = — 

\ ^ , a3 l 

1(1 -f) 2 J 

r 

er 


Here p' represents the momentum in the primed coordinate system with 0, IM H3 

= ^ 


( 8 ) 

(9) 


and 


a e A s = — 
er 


±>i 

(1-i) 


i<j2 
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(10) 
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Using Eq.© and Eq. dlOl) . the Dirac equation in the presence of the conical defect can be written as 


[H — M]A> = 0, where H = a 3 
From Eq. m we get 




H 2 - M ' 2 = id 2 + -a r - (. M 2 - u 2 ) + 


2Ea 


where 


^ = _H^±f) 2 


( 1 -t ) 2 (i-l) 


(-ids ± f) 1 . 

V3— - ^r~ + {a -p-zaa 2 . 


The Green’s function in the presence of the conical defect can thus be written as 

p oo ^ r~— 

G(r, fo; E) = —i[H + M] / dse irs{li2 - M2} J — 5{r - r 0 )S(d - 6 0 ). 


Let | Z(6) > be the normalized eigenfunction of the operator Z. 

Z\Z K (6) >= —k 2 \Z Kj {6) >, 


( 12 ) 


(13) 


(14) 


(15) 


where — n 2 denotes the eigenvalue of Z. \Z(9) > has two components corresponding to the two sublattices of the 
graphene. Now consider the projectors [lj|[l8|] D K (0,0 O ) = 1-2 K (0) >< Z K (6 0 )\ which satisfy the relations 


and 


ZD K (Mo) = -K 2 tl K {0,0 o ) 


Y / n K (9,0 o ) = 6(0-e o ). 


(16) 

(17) 


The eigenvalue — k 2 is determined from Eg. (1161) . which can be written as 


( -id e ±^) 2 _ (-id e ±$) o _ 1 , 2 

(l-#) 2 (1-f) +Q! * +h 

a 


—a 

(-i0e±t) 2 I (-i9s±f) , 2 1 , 2 

' + (1-f) +« -4+ K 


^12 

^21 ^22 


= 0. 


(18) 


Using the completeness of the total angular momentum eigenfunctions, the components of the projectors can be 
expressed as a linear combination of the functions e L ^ e ~ e °^ and e~ l ^ e ~ s °\ where j is taken to be a positive half 
integer. Using this and solving Ea. (fT51) we finally get k = (77 ± |) with ?y = \Jv 2 — a 2 , v = and 


n 


n 


= (v + v) c ij(g-e 0 ) + iv ~ v) c -ij(e-e 0 ) = 


47177 


Amj 


and ill, = - [ e ij(»-6o) +e -me-B 0 )-i = q 

Ami 


(19) 


Note that these projectors are different from what is obtained for planar graphene. This is due to the presence of the 
quantity v in the projectors which explicitly depend on the number of sectors n removed from the graphene plane to 
introduce the conical defect. The Green’s function in the presence of the topological defect can then be written as 


( 20 ) 


G{r, 7"o; E) = -i a 3 (e + ^ ) (~id e ± %) ~ *i ^ j + -A 

x VfT(Mo) / 00 rf S e 2i8Ba - 2i ^K-^ + r]+(^ 2 ^ 2 )i} & (r _ r o). 

„ Jo V r 0 


The Green’s function explicitly depends on the conical defect through the appearance of the quantity n in it and it 
is different from the Green’s function for the planar case. 
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FIG. 2: Dependence of polarization charge on subcritical Coulomb potential using Green’s function. 


In order to proceed we follow the operator technique for Green’s function in QED developed in [18| . which was 
adapted to the problem of planar graphene in [20|. The key point in this approach is to note that the operators 


-rdl - d r + ^ 


, r and —* ( rd r + i) form the generators of the 0(2,1) algebra . Using the commutation relations 

of these 0(2,1) generators, the action of the exponential in Ea. (l20l) on S(r — ro) can be evaluated. The integral 
representation of the Green’s function can be written as 


G(f= i 


03 




tO 2 




r{ 1 -t) 

I< 


( K -id e ±l)-o 1 U + ±) + 


ds 


o sin (Ks) 


M 


exp [2isEa + iK(r + ro) cot(A's) — inn] J^k, 


( 2Ky/rro\ 
\ sin (Ks)) 


( 21 ) 


where K = ^J(A4 2 — E 2 ). The Green’s function in Eq. (l2lT) has been expressed in terms of certain operators and 
an integral. The information about the conical defect appears in the Green’s function through n, n and G K . When 
n = 0, Ea. (l2Tl) reduces to the Green’s function for planar graphene. The crucial point to note at this stage is that 
the integral in Eq. m lias the same analytical structure as that for the planar graphene, which is a consequence of 
our treatment of the conical defect. Due to this, we can follow the analysis in [20] in order to evaluate the integral 
in Eq. (l2lT) . This calculation involves the choice of a suitable contour of integration [l9|, [20]. The resulting expression 
for the induced charge is divergent, which requires a renormalization as in the case of QED |l9]. The condition that 
the total induced charge is zero is imposed as a physical requirement in the renormalization process. This is possible 
due to the presence of the infrared cutoff given by the quasiparticle mass M and it is for this reason that M was 
introduced. The details of the calculation in the planar case can be found in [19, 20], following which, we get the 
induced charge density in the M —> 0 limit as 


Pi(r) = QG$(r), 


where the induced charge Qg is given by 


Qg = -e 


7rcr 8 V—' T I tLX 

-b - > Im { — 

9 IT ^ 


ia ln (?7 — ia) 


2v 


+ lnT(7y — ia) — (rj — ia)t/j(ri — ia) — iaz't^z') 


( 22 ) 


(23) 


where —e denotes the electronic charge and 0 and tp' denote the digamma function and its derivative respectively. 
Thus, for a subcritical Coulomb charge in gapless graphene with a conical defect, the polarization charge is localized 
at the apex of the graphene cone, which is also the position of the external Coulomb charge. 

In order to illustrate the effects of the conical defect, in Fig.© we have plotted the dependence of the polarization 
charge \Qg\ (in units of electronic charge) on the Coulomb charge strength a. The plot shows the behaviour of the 
induced charge for both n = 0 (planar graphene) and n = 1, which corresponds to the conical defect. For n = 0, the 
plot shows a marked deviation from linear profile at around a = 0.5, which corresponds to the critical value of the 
charge in the planar case. For n = 1, the plot shows a similar deviation starting with a smaller value of a = 0.3. This 
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behaviour is consistent with the fact that the critical charge in the presence of the conical defect depends on n and 
that for n = 1. the critical charge has a value given by 0.3 [38|. To get an analytical feeling for this phenomenon, it 
is useful to consider the expansion of Qg as a function of a. For n = 1, we get 


Qg = -e 


— + 5.289a 3 + 27.437a 5 


(24) 


For the planar case, the leading term has the same value but the coefficients of a 3 and a 5 have values 0.783 and 1.398 
respectively j2(| . From these values it is clear that the deviation from linearity sets in at a lower value of a in the 
presence of the conical defect, which is a clear indication of the influence of the topological defect on the magnitude 
of the induced charge. 


B. Induced charge from Friedel Sum Rule 

The Friedel sum rule for a gapless graphene cone is given by 

A N = -J2[SAEf)}, (25) 

7T z ' 

3 

where AN represents the change in the number of charged particles due to the Coulomb potential, Sj represents the 
scattering phase shift in the j-th angular momentum channel and Ep is the Fermi energy [27]. For planar gapless 



a 


(b) 



FIG. 3: (a)Dependence of on subcritical Coulomb potential is shown for j = \ and § when n = 0 and 1. (b)Dependence 
of polarization charge on subcritical Coulomb potential is shown using both Friedel sum rule and Green’s function technique. 
We have considered the effect of j — \ only. The solid lines correspond to Green’s function technique and the dashed lines 
correspond to Friedel sum rule. In this paper the positive values of a correspond to an attractive Coulomb potential. 


graphene, the induced charge from the Friedel sum rule was discussed in Q in detail. Here we will discuss the problem 
in the presence of the conical defect. The expression of scattering phase shift for the massless graphene with a conical 
defect in presence of a subcritical Coulomb charge is given by 

Sj(k) = asgn(E) ln(2fcr) - sgn (E) ^argf ^(v - ia)\ + arg[r(l + rj + ia)]^) - ^(77 - \u\), (26) 

where k = —E. The positive energy values correspond to the electrons whereas negative energy values correspond to 
the holes in our work. The detail of the calculation can be found in Appendix A. From Ea. (E^l) we can see that the 
scattering phase shift contains the term v which depends on n. Note that the scattering phase shift is asymmetric 
with respect to the sign of a. For planar graphene, the phase shift has been discussed in [3, [f§, [48[ in detail and we 
recover the same result when n = 0 . 

The polarization charge is now defined as Qf = —eAN. We have shown the effect of two different angular 
momentum channels j = ^ and § on the scattering phase shift in Fig.([3Ji). The plot shows that only in the lowest 
angular momentum channel j = i the sign of 8j has an expected nature for attractive and repulsive Coulomb potential. 
This feature has already been noted for planar graphene | 8 ] and continues to hold in the presence of the conical defect 























7 


as well. Due to this reason, we consider the effect of the lowest angular momentum channel only when calculating 
\Qf\- In our plot we have neglected the customary asgn (E) ln(2 kr), which is due to the Coulomb potential 3]. 

In Fig. (|3 Jd) a comparison between \Qg\ and |Qf| is shown. For n = 0, both of them have similar behaviour for 
low values of a. In the presence of the conical defect, the difference is apparent from a « 0. It is known that there 
are subtleties associated with the Friedel sum rule in the presence of singularities [39f. Whether similar effects are 
playing a role in the presence of a topological defect requires further investigation. 

We would now like to study how the conical topology affects the physical observables such as the local density of 
states (LDOS) and the conductivity in graphene. The LDOS /x is defined as 

LL=-J2\*(E,r)\ 2 , (27) 

7r L ' 

3 

where 'k(£ , ,r) represents the normalized eigenstates directly obtained from the solution of the Dirac equation. We 
have plotted the dependence of fj, on the energy of the quasiparticles in a graphene cone in Fig.(@]r) where the distance 
from the Coulomb charge has been chosen as r = Iq where Iq is of the order of the lattice scale in graphene and 
energy is given in units of Iq . When a positive Coulomb charge is introduced in our system it attracts negative 
charge carriers and repels the holes. Asa result, the density of negative charge carriers increases and that of the 
positive charge carriers decreases near the impurity. The same effect can be observed from this plot. In the presence 
of a conical defect with n = 1, the plot shows that the density of both positive and negative charge carriers decrease 
compared to the planar case. It may be possible to detect this effect of the conical defect with the STM. 



(b) 



EPo -1 ) a 

FIG. 4: (a)Dependence of LDOS on energy of the quasiparticles in the subcritical region. We have considered r = lo in our 
plot. (b)Dependence of or on Coulomb potential strength is shown. 


The transport cross section cr^ rans is related to the scattering phase shift by the relation 


CT trans ~ j. X/ 


sin 2 {Sj - <5j+i) 


4 

r T> 


( 28 ) 


where k = —E and Sj is given by Eq . (l26D for the subcritical regime. The conductivity is inversely proportional to 
the transport cross-section SLUSHa l48j |. We have plotted the dependence of ctt on subcritical Coulomb potential 
strength in Fig. (Or). The plot is asymmetric with respect to the sign of the potential because the scattering cross 
section of conduction electrons depends on the polarity of the external Coulomb charge. We can see from the plot 
that when the value of n increases from 0 to 1, the value of <jt also increases. This can be interpreted to be arising 
from additional interactions associated with the conical defect, which leads to a decrease of the conductivity. 


III. SCREENING OF SUPERCRITICAL COULOMB CHARGE IN GRAPHENE WITH CONICAL 

DEFECT 

For planar graphene, the induced charge density in the supercritical regime has been calculated using both scattering 
phase shifts Q as well as the Green’s function [4Ci(. Here we analyze the effect of conical topology on the induced 
charge in the supercritical regime using the scattering phase shift. In the presence of a conical defect, the supercritical 




















regime sets in when a > v for given values of n and j [38| . Defining /3 = ir\ = \Ja 2 — v 2 and using zigzag edge 
boundary condition, the scattering phase shift can be written as [38j 


where k = —E, 


Sj(k) = arg[e is S n(E)4(fc) + ce - ^ 11 ^)^)] - sgn(£7)arg(/i ai p) + asgn(F) ln(2£r), 

c = h a>p = C = 715 and e2i?(fc) = ^§e 2 ^n(- 2 sgn (E )« 0 ). 


(29) 

The 



FIG. 5: Using Friedel sum rule dependence of polarization charge on wavenumber klo is shown in the units of electronic charge. 
Here we have considered the effect of the holes (sgn(I?) = —1) for the lowest angular momentum channel j = \ where a = 1. 

scattering phase shift in the supercritical region has a strong energy dependence through the first term in the right hand 
side of Eq. (EH1) . As discussed in [3], the induced charge within a radius r is not affected by modes with wavelengths 
larger than r. Thus we get 


—4p 1 

QM = — 

V<OL 


(30) 


Neglecting the customary asgn(E?) ln(2/cr) term of Eq. (l29l) and using Eqns. (l29l) and (IdOl) . Q(r) can be expressed in 
the units of electronic charge as 


Q(r) = 


4sgn(A) 


E 


is<ct _ 


/3 ln (irr)+tan 1 


_ Fi{a, v) cos[/(a, v) - 2sgn(£)/31n ((&)]_ 

F 2 (a, v) + sgn(E)Fi(a, v) sin[/(a, v) - 2sgn(E)/3ln ( l f)} 


+ F 3 (a,is) 


, ( 31 ) 


where F\(a,v) = e s g n ^) ,r/3 C|r(l + i/3 — ia )|, F 2 (a,v) = |T(1 — i/3 — ia )|, F 3 /a,v) = —tan -1 £ + arg (/i ajj g) and 
/(a, v) = —2sgn (E) tan -1 ( + 2argT(l — 2i/3) + argT(l + i(3 — ia) — argT(l — */3 — ia). In Fig.©, using the Friedel 
sum rule, the dependence of |Q(r)| on the wavenumber klo is shown in the units of electronic charge. We consider the 
effect of the holes in our plot for which sgn(A) = —1. The kinks appearing in the plot correspond to the formation 
of quasibound states. 

From Q(r) we can calculate the induced charge density p/(r) in the supercritical region of the graphene cone using 
the relation [s] 


Piir) = 


1 dQir) 


2irr dr r 2 

The coefficient of the power law tail of induced charge density is given by 


A,(r) = % 

7 T z 


sgn(£)-2E 

^2 


(32) 


(33) 


where 


n l ^ 2 
D i = I — sec 


/ + 2sgn(F)/?ln ( — 


tan 


/ + 2sgn(F)/31n ( — 


sgn(F) sec 2 / + 2sgn(A)/31n ( — 
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and 


D 2 




f + 2sgn(E)/3ln 



+ sgn(E) tan 


/ + 2sgn(£’)/31n 




The summation is restricted to those values of v such that a > u, which is needed for supercriticality. From Eg . (1531) 
we can see that D„(r) is a log-periodic function of the distance from the Coulomb impurity. The period of variation 
and the average value of D v {r) are given by and — 2S g I F^ £ )' 3 respectively. These two quantities depend on the 
sample topology through the term /?. For n = 0, the period of variation of D v (r) agrees with the result obtained for 
planar graphene using the Green’s function technique j4(j and the average value of D v (r) matches with the constant 
coefficient of the power law tail given in [3]. 



e(/<t‘) 


(b) 



FIG. 6: (a)Dependence of standing wave oscillations in LDOS on energy of the quasiparticles in the supercritical region. We 
have considered r = 101o in our plot. (b)Dependence of or on Coulomb potential strength. We have considered klo = 10 -5 in 
the plot. 


The dependence of the LDOS n on the energy of the quasiparticles has been plotted in Fig.([(3i). The LDOS has 
been calculated using eigenstates obtained from the solution of the Dirac equation [38| , using the angular momentum 
channels for which a is supercritical. The amplitude of the oscillations in LDOS depend on the topological defect and 
they decrease for n = 1 compared to the planar case. 

The transport cross-section of graphene with a conical defect has been calculated using the relation given in Eq. (E51) , 
where the scattering phase shift Sj in the supercritical regime is given by Eq. (1291) . The dependence of <jt on external 
Coulomb potential strength is shown in Fig. ([HJo) •> for both sub and supercritical regime. The peaks in this plot 
correspond to the quasibound states which occur when the charge is supercritical. It is seen that the peaks start 
appearing for a lower value of a when n = 1 compared to when n = 0. This is consistent with the fact that the critical 
charge for n = 1 has a value 0.3 compared to the value of 0.5 for the planar case. As in the subcritical case, the plot 
for <jt is asymmetric with respect to the sign of the potential as the transport cross section of conduction electrons 
depends on the sign of the external Coulomb charge. 


IV. GENERALIZED BOUNDARY CONDITIONS 

The quasiparticles of graphene obey Dirac equation in the long wavelength approximation. On the other hand, 
the external Coulomb charge and the topological defect can interact with the graphene lattice and give rise to short 
range interactions. Such short range interactions cannot be directly included in the Dirac equation, which is valid 
only in the long wavelength limit. One way to include the average effect of such short range interactions is through 
the choice of appropriate boundary conditions [43!]. A technique to classify all such boundary conditions is provided 
by the method of self-adjoint extensions due to von Neumann [43], 0, 0, (47| . This formalism provides all allowed 
boundary conditions which are consistent with probability current conservation and unitary time evolution. 

The generalized boundary conditions are relevant for graphene with an external charge and a conical defect only 
when 0 < [45]]. In this case, the allowed boundary conditions are labelled by a single real parameter ui defined 

mod 27 t, also known as the self-adjoint extension parameter. The precise value of this parameter cannot be predicted 
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by theory alone and has to be determined empirically. For this class of generalized boundary conditions, the scattering 


(a) 



(b) 



a a 

FIG. 7: (a)Dependence of polarization charge on subcritical Coulomb potential using Friedel sum rule and generalized boundary 
conditions is shown. The solid lines correspond to cj = 0 and the dashed lines correspond to w = ^. (b)Dependence of or on 
Coulomb potential strength is shown in presence of generalized boundary conditions. The solid lines correspond to cj = 0 and 
the dashed lines correspond to uj = ^. 


matrix S and the corresponding phase shift 5{k) are given by 


s = e 2iSgn(E)5(fc) 


J/e 2iSgn(E)aln(2fcr) 


A + B 
v 2 A* A B* 


(34) 


where 


A = 


*7rSgn(B)r, 

e 2 


r(i + t] 


( k) 2ri e 


-i^sgn (E) V 
2 


ia)r(l — rj — ia) 


(35) 


iirSgn(E)„ 

e 2 (—rj + ia) 


F(1 + r] — ia)r(l — rj + ia) 


(k) 2r, e 

r(i - v 


in sgn t E)„ . 

2 (ry + ia) 


fa)r(l + 77 + ia) 


e 2 . 


(36) 


See Appendix B for the details. It is clear from Ea. flMl) that the scattering phase shift depends on the self-adjoint 
extension parameter oj explicitly. 

The induced charge can be calculated from the phase shift using the Friedel sum rule. In Fig. ([TJa,) we have plotted 
the induced charge as a function of the Coulomb potential strength a with the condition that 0 < 77 < A. The 
polarization charge depends on the conical defect as well as on the value of the self-adjoint extension parameter ui. 
The non-zero value of the induced charge for n = 1 and a = 0 is indicative of the additional interactions due to 
the topological defect whose average effect is captured by the generalized boundary conditions. In Fig.Qj) we have 
also plotted the dependence of ctt on a and ui for n = 0 and 1. The contribution of only j = y channel is relevant 
when 0 < 7j < -r. It is seen that the effect of the generalized boundary condition is more pronounced when n = 1 as 
compared to the planar case. This is consistent with the fact that generalized boundary conditions are more likely to 
be important in systems with topological defects. 


V. SUMMARY 

In this paper we have analyzed how a conical defect in gapless graphene affects the Coulomb screening and related 
physical observables. In the subcritical regime, we have used an operator approach to the Green’s function in order 
to evaluate the induced charge. Our calculation shows that the Green’s function (l2lT) in the presence of the conical 
defect can be expressed as a product of certain operators and an integral, both of which depend on the conical defect 
through their dependence on n. The evaluation of the integral in (1211) requires a careful choice of the contour and 
a proper renormalization procedure. In our way of treating the conical defect, the integral in (1211) has the same 
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analytical structure as that for the planar case, although they are not identical. This observation allows us to write 
the result of the integral in (l 2 lT) from the knowledge of the corresponding integral in the planar case, leading to a great 
simplification of the problem. The final result for the induced charge Qq in (12311 explicitly depends on n through 77 . 
In other words the induced charge depends on the number of sectors that have been removed from a planar graphene 
lattice in order to form the graphene cone, or equivalently it depends on the sample topology. The n = 0 limit of (l23l) 
agrees with the result obtained for planar graphene [ 20 ]. 

Another interesting effect of the conical topology can be seen from Fig. 2, which shows the dependence of the 

induced charge on the subcritical Coulomb strength a. It is clear that for n = 1, the nonlinearity sets in much before 

that that for n = 0, the latter being the planar case. This onset of nonlinearity happens as the Coulomb strength a 

approaches the critical value. This plot therefore shows that the critical value of the external charge actually depends 

on the topological defect. This striking interplay between dynamics and topology can be understood by noting that 

unlike the planar case, here the critical charge is governed not just by the angular momentum j but by the quantity 

v = {HJ. The behaviour of the induced charge, which is a physical observable, provides a further evidence of 
6 

this interplay. 

In the subcritical regime, we have also calculated the scattering phase shifts and have obtained the induced charge 
using the Friedel sum rule. For the planar case, these two approaches give similar values of the induced charge for 
small values of the external Coulomb charge and the results begin to diverge as the critical value is approached. In 
the presence of a topological defect, this divergence sets in for a much lower value of a. This is possibly related to the 
subtleties associated with the Friedel sum rule in the presence of singularities [39!], although the exact reason behind 
this observation requires further investigation. For a given value of the subcritical charge, the transport cross-section 
has a higher value in the presence of the conical defect. Similarly, the LDOS in the subcritical regime also reflects the 
effect of the conical defect. 

In the supercritical regime, again the induced charge as well as the LDOS and the transport cross-section depend 
explicitly on the conical defect. For the planar case using a Green’s function approach, it has been observed that 
the polarization charge density in the supercritical regime exhibits a power law tail. The coefficient of the power law 
tail can be written as a sum of functions which vary log periodically with the distance from the Coulomb impurity 
[40l | . We have obtained the induced charge in this regime using the scattering phase shift and the Friedel sum rule. 
What we find is a similar power law tail where the period of the oscillation depends on the sample topology through 
its dependence on n. Moreover, the planar limit of this quantity as obtained from our calculation agrees with that 
obtained from the Green’s function approach [Zo[. The transport cross-section in the supercritical regime shows that 
the Fano resonances UM begin to appear for a lower value of the Coulomb charge in the presence of the conical 
defect. This is consistent with the observation that the critical charge is modified by the conical defect. 

Finally, for certain ranges of the system parameters the graphene cone with an external Coulomb charge admits 
generalized boundary conditions mum. These boundary conditions incorporate the effect of any short range inter¬ 
actions that may have been caused due to the interaction of the external charge as well as the topological defect with 
the graphene lattice. Such interactions cannot be directly incorporated in the Dirac equation which is valid only in 
the long wavelength limit. One possible way to take into account such effects is through these generalized boundary 
conditions. We have shown that the observables such as LDOS and the transport cross-section depends both on the 
sample topology as well as the choice of the generalized boundary conditions. 


Appendix A: Derivation of scattering phase shift for subcritical region in graphene cone 


In this appendix the expression of scattering phase shift for the massless graphene with a conical 
of a subcritical Coulomb charge is obtained. The Dirac equation for the massless graphene cone is 


defect in presence 
given by [381 ] 


where 


We assume 


H 'F = E 'F, 


H = 


— d r — 


rC-f) 


dg ± 


A 1 -#) 


2 r 



a 

r 




g-*£> r 7j-(l/2) e ij8 


(Al) 


(A2) 
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For doing the calculations we define two functions U^\r) and V^\r) by 

^ , A’ > ( r ) = V^{r) +U^(r) and ^^(r) = V^(r) — U^\r) 

. Using these functions from Eq. m we get 

cPV^tx) dV^\x) . . . 

' + (1 + 2 ?y — x )—--(ry + ta) V' 3 '{x) = 0 , 


da ; 2 


da; 


where x = —2ikr, with k = —E. 

The solution of Eq. (IA3I) which obeys regularity at the origin is given by 

V U) {x) = CM (ry + ia, 1 + 27y, x ). 


(A3) 


(A4) 


Here M(a, b , a;) is the confluent hypergeometric function [49}] and C is a constant which depends on the energy of the 
system. 

From Ea. (IA4l) and Ea. dAll) we have 


lf{j) (a;) = C ^ ^ M (1 + ?y + ia, 1 + 2?y, a;). 


(AS) 


Using the asymptotic form of (a;) and U^\x) we can find out the expression of the scattering phase shift to be 


Sj(k) = -sgn (E) 


—aln(2 kr) + arg[r(l + 7y + ia)] — — tan 


-l 


( 4)1 

Trry 

IS7T 

\vj \ 

2 

2 


Ea. dA 6 H is same as Eo. (l26l) . For n = 0 this result reduces to the one given in 48). 


(A 6 ) 


Appendix B: Derivation of scattering phase shift with generalized boundary conditions 

The radial part of the Dirac operator H in Eo. dA2l) can be written as 

Hrad — 

The equation which determines the domain of self-adjointness of H ra d is given by 


-{<9 r - {u- 4)4} 


R+ (*+!)£} 


H ld xl, ± = ±J*±, 

where 1 is a constant with dimension of length. It can be set equal to unity with a choice of suitable units. 




V£ i} (r) + U^(r) 


U), 


i(V±\r) -U'i’ir)) 


(j), 


t jx 2 _ 


(Bl) 


(B2) 


(B3) 


The total number of square integrable, linearly independent solutions of Eauation dB2l) gives the deficiency indices n± 
phi ]. which gives a measure of the deviation of the operator H ra d from self-adjointness. When n + = = 0, H ra d 

is essentially self-adjoint in T>o(H ra d) where V o = Cfi°( 1Z + ) consisting of infinitely differentiable functions of compact 
support in the real half line 1Z + . When n+ = n_ ^ 0, H ra d is not self-adjoint in T>o(H ra d) but it can admit self-adjoint 
extensions. When n .|_ 7 ^ n_ H ra d cannot have self-adjoint extensions [46} . 

The solution for 4^^+ is given by 


'&A+ = e t 
As r —> 00 , 4>a+ —> 0 and when r —> 0, 


2r\ 1 / 2 r 

U ( 1 + 7 y - ia, 1 + 277 , — J - U lij — ia, 1 + 27?, — 


f 7 ! 2 


J \'&A+\ 2 'rdr ~ J r 2r, dr + converging terms. 


(B4) 


(B5) 
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Thus \I/a+ is a square integrable function for the range 0 < 77 < Similarly the entire radial wave function can 
be shown to be square integrable for the specified range of 77 and the deficiency index n + = 1. Similarly we can 
show that n_ = 1 when 0 < 77 < Thus for 0 < 77 < 4 we get n+ = n- = 1 and the radial Hamiltonian H ra d 
admits a one parameter family of self-adjoint extension 46]. The domain of self-adjointness of H ra d is given by 
T^ui{H ra d) = 'D Q ( y H ra d) ® {e 1 ^ 4/+ + where w £ R mod 2n is the self-adjoint extension parameter. 

As r —> 0, 


^A + = 


j/sin7r(l + 2 77) 


(v + 77 + ia) 


and 


tE 'A. = 


T (1 — 77 — 7 o)T(l + 277) 
(v + 77 + ia) 


2 — — 


-2 jj 


(v — 77 + ia) 


I 7 sin 7 r(l + 2rj) 


T (—77 + 7o)T(l + 2 77 ) 


r 2 — - 


T(1 + t? - *o)T(l - 277) 

“ 2rj O' - V +H 

r( 7 ? + *a)r(i — 277 ) r 


(B 6 ) 


(B7) 


Note that the physical scattering states are given by the solution of Eq. (ED, from which we get 


n-bjkrcn 0 / + * a ) 


-2n ( _7 7 + *«) 


47,4(0;) = r v 2 e' ,Kr {C\X-!—L --M (1 + 77 + ia, l + 2rj, x)+C2X 2r] - — - -M (1 — 77 + ia, 1 — 277, x) 

+ CiM(rj + ia, I + 277 , *) + C 2 X~ 2v M (—77 + 70 , 1 — 277 , a;)}. (B 8 ) 

We now match the behaviour of the physical wave function with a typical element of T> u (H ra d) as r —> 0. In the limit 
r —> 0 Ea. dBSD gives 

(v — 77 + ia) 


T / x „ 72 + 77 + 70 1 

T A (r)=Ci--- r" 2 +C 2 


-(—2 ik) 




77 77 

In the same limit, the upper component of a typical element of the domain T> ul (H ra d) is given by 

lim A(e^47,4 + + e^~^>A-)- 

r—> 0 

Comparing m and (iBTOl) and using (1B6I) and dB7l) we get 


Ci = 


A7T 


and 


sin7r(l + 2rj) 

C 2 = -(- M) 2r > 


e 2 


e 2 


T(1 - 77 - 7a)T(l + 277) r(-77 + 7a)T(l + 277 ) 


Xn 


sin7r(l + 277) 


e 2 


+ 


e 2 


When r -A- 00 we note that 

vMO = (-2 ik)-n(-iy 
+ (-2 iky^-iy 


T(1 + 77 - *a)T(l - 277) r(r? + *a)r(l - 277) 

g— i[kr— a ln(2fcr)] 


(B9) 

(BIO) 

(Bll) 

(B12) 


rj + ia T (1 + 2rf) -77 + ia T (1 — 2if) 

Ci- , — , . , +C 2 


Ci 


77 r(i +77 + 70;) 

T(1 + 277) 


77 r(i — 77 + 70) 

r(i-2i) 

T (1 — 77 — 7 a) 




e ^7^(r?+^Q:) | ^ V ( — itv(— r)+ia) 


a i[kr—a. ln(2fcr)] 


V~r 


(B13) 


T(1 + 77 — ia) 

Substituting the expressions of Ci and C 2 from Ea. dBllI) and Ea. dB12l) in Eci. dB13l) the scattering matrix S and the 
corresponding phase shift S(k) are obtained as 

A + B 


S = e 27Sgn(B)5(fc) = _^ e 2iSgn(B)aln(2fcr) 


where 


A = 


■ *»sgn (J ), 
e 2 — (kl) v e 


1 A A* + B* 

-7„sgn(E)„ \ 


T (1 + 77 — 7 o)T (1 — 77 — ia) 


(B14) 


(B15) 


B = 


ixr SgIl(E)„ 

e 2 (—77 + 70 ) 

T (1 + 77 — 7o)T(l — 77 + 70 ) 


(fcQ^e ~" Sg 2 n(E),, (77 + 7o) | 
T (1 — 77 — 7 o)T(l + 77 + ia) j 6 


(B16) 


It is clear from Eq. dBldj) the scattering matrix and the phase shift explicitly depend on the self-adjoint extension 
parameter oj. Ea. dBllI) . dB15l) and dB16l) are same as Eg. (1+0) . (l35l) and (I5S1) given in section 4, with the choice of the 
dimensionful parameter 1 = 1 . 
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